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Abstract 



(N 

^ijT . The field-theoretical description of quantum fluctuations on the 

Q"^ I background of a tunneling field a is revisited in the case of a functional 

Schrodinger approach. We apply this method in the case when quan- 
tum fluctuations are coupled to the a field through a mass-squared 
Qh! term, which is 'time-dependent' since we include the dynamics of a . 

(— I I The resulting mode functions of the fluctuation field, which determine 

the quantum state after tunneling, display a ,previously unseen, reso- 
nance effect when their mode number is comparable to the curvature 
kS I scale of the bubble. A detailed analysis of the relation between the 

d ' excitations of the field about the true vacuum (interpreted as particle 

creation) and the phase shift coming from tunneling is presented. 



1. Introduction. 

False vacuum decay through barrier penetration has been studied extensively 
by a number of authors [1-10]. However there is a lot to be understood yet 
about tunneling phenomena, especially if one considers its importance as 
related to the very early universe. Tunneling is often related to the existence 
of a potential barrier ,but while in quantum mechanics (QM) that barrier is 
in real space, in quantum field theory (QFT) one has penetration from one 
field configuration to another. In QFT the effect of a tunneling potential 
can even be produced through the spatial gradient of the field. Thus the 
potential barriers in (QM) and (QFT) have a different nature. 

It is well known that tunneling in QFT is a large fluctuation or instability 
of vacuum, a classically forbidden process where energies are not eigenstates 
of the hamiltonian due to barrier penetration. Since particle creation is 
another phenomena related to vacuum fluctuation [11-18] it is of interest 
to probe tunneling via studying the relation between particle creation and 
tunneling. 

The tunneling amplitude depends strongly on h since it is the exponen- 
tial of the negative Euclidean action divided by h, while particle creation 
is insensitive to it. For more details see [7,9]. Thus, the question we ask 
here is: What is the relation between tunneling and particle creation in the 
dynamic background of the bubble that forms in false vacuum decay? We 
will consider the functional Schrodinger equation for a quantum field $ in a 
tilted quartic potential V{^) (fig. 1). The field is split into a classical part, a 
(hereafter the tunneling field) undergoing phase transition from false to true 
vacuum, and the fluctuation field, ip around a, giving rise to particle creation 
through coupling to a with a mass-squared term, 'm^(o").The classical field 
a is defined as usual to be the one that extremizes the action. In this case, 
assuming an WKB ansatz for the wavefunctional, (Section 2 reviews this for- 
malism), is constistent with the definition of a, since the expectation value 
for the lowest order WKB type solution recovers the classical solution a in 
the appropriate limit. However generally this statement would not be true 
particularly in the cases when there are more than one degree of freedom and 
the WKB branches are not decohered. The multidimensional tunneling in 
WKB formalism was initially developed by Banks, Bender and Wu [5] and 
later revisited and improved by T. Vachaspati and A. Vilenkin [6,7] as well 
as the Kyoto group [Hamazaki et al 8,9]. 



The reason for using this method is twofold: 

- It allows us to relate the dominant escape path,DEP [22], of a to the the 
usual formalism of particle creation [11,12], i.e. relate the fluctuations ip 
around the DEP of the o"-field to the particles produced as the result of tun- 
neling. 

- While the approach is similar to QM (i.e. solving a Schrodinger equa- 
tion) there is a clear relation with the second quantization picture (shown 
rigorously by [2,3,7,8]). 

The method follows closely [7,9] and we review the formalism in Sect. 
2. The ansatz of the wavefunctional is assumed to be a superposition of 
Gaussian packets centered around the classical vacuum outside the barrier 
and around the DEP under the barrier. One then solves the functional 
Schrodinger equation for the a-field and the fluctuation (/9-field on the back- 
ground of o", with the appropriate boundary conditions at the inner turn- 
ing point. The wavefunctional that describes the quantum state of the two 
fields, i.e. that gives the amplitude for the fields to be in configuration a, 
ip, is a 1-parameter family of solutions (the parameter can be thought of 
as a time-coordinate) and is analytically continued from under the barrier 
through the turning point, to the region of the true vacuum. The solution for 
the fluctuation field is given in Section 3. When the solution is analytically 
continued, one ends up with a squeezed state of mixed positive and negative 
frequencies. The coefficient of the negative frequency component than allows 
one to determine the particle creation number in each mode. 

We find an interesting new feature for the fluctuations, namely resonance, 
by taking into account the dynamic of the background fleld a. A discussion 
of the results and open issues to be tackled in the future is given in Section 
4. 

All along we have assumed the 0(4)-symmetry and the thin wall approx- 
imation. In the appropriate limit our result recovers those obtained by the 
Kyoto group [8,9]. 

2. Review of the Functional Schrodinger FormaHsm 

The application of this formalism in connection with vacuum tunneling in 
QFT was initially developed by Banks, Bender and Wu [5] and later on 
revisited and improved by [2,3,7,8,9]. We consider a scalar field in Minkowski 



spacetime with the action 



S= (fx 



where 



■{d.-bf - V(») 



n*) = 5(*^-<.T-^(* + a). (^«1 



[see Fig.l] 



Fig.l The Potential for False Vacuum Decay via Tunneling 
We want to study the fluctuation field around the dominant escape path 
(i.e. around the tunneling field a). Thus we split the field into the classi- 
cal/tunneling field "a" and fluctuations around a, denoted by y?, i.e. $ = 
a + (p . 

Under the barrier, a interpolates between the 2 vacuum states by taking 
the dominant escape path. The solution for the a-field in the case of potential 
V{a) is studied by many authors through the use of instantons in the dilute 
gas approximation under the barrier, or by means of WKB and loop expan- 
sion methods outside the barrier.lt has been shown that both perturbation 
methods (dilute gas of instantons inside and 1-loop approximation outside 
the barrier) match consistently at the turning point, only when the continuous 
symmetries of the problem (in our case 0(4)) are the same under the barrier 
and outside, where one is finding solutions around a ~ cr±. outside) [2,3]. 
This matching procedure can be quite complicated if one breaks different 
subgroups of the continuous symmetries inside and outside the barrier. 

2.1 Particle Creation as excitations of the tunneling field. 

Consider a scalar field theory $ in a quartic potential 



As noted earlier, we split the fields in 2 parts: a, the classical field 
( where ^|$=o- = 0) and ip,the fiuctuation field. The false vacuum is 
located at (T_ ^ — a— g^ + 0(e^), where V{a_) ~ + O(e^). The true 
vacuum is located at 0-+ ~ a — g^ + O(e^) , where V{a^) ~ — e + O(e^) 
. Then the Lagrangian for the (a + (/9)-fields becomes 

L = [lid.af - V{a)] + i {av^)2 - m\a)v'} = L^ + L, (2.1) 

where m'^{a) = ^^ and (□(/? — m'^(a)(p) = We have ignored higher 
order terms in the potential that would have contributions of order ip^ 
and higher. This is a consistency requirement for matching the two 
different perturbation methods at the turning point. When we solve 
the functional Schrodinger equation, the solution for the a-field under 
the barrier will contain terms up to {h~'^) coming from the zeros in 

the determinant factor K = (^j det""*^'^ | — (9^(9^ + V"{a)\, where 
5*0 is the action of the bounce in 4-dim. These zeros are related to 
the 4 translational symmetries of the problem, and it has been shown 
(see Coleman [19]) in the dilute gas approximation that each symmetry 
brings a {h~ ' ) factor, so counting zeros in K is equivalent to counting 
powers in the coupling constants. Further, Gervais, Vegas, Sakita, 
[2,3] have shown that for the problem at hand, expansion around the 
vacuum states, outside the barrier, which in our case corresponds to 0"+, 
matched (analytically continued) with the solution under the barrier 
is in complete agreement with the multi-instanton (dilute gas approx.) 
approach [1]. Thus we can keep terms up to (p"^ in the Lagrangian and 
ignore higher order terms. 

We work with the functional Schrodinger equation and use the results 
of Vachaspati and Vilenkin [7] for the background field a. 

2.2 Functional Schrodinger Equation and Solution for the Classical Field. 

From the Lagrangian in (1) we can write 

H^H^ + H, = [\nl + V{a)] + IttJ + \m\a)^\ 

The quantum state of the 2 fields is described by a wavefunctional 
■?/'[(t(x, t), y9(x, t)] which gives the amplitude for the fields to be in con- 



figuration o"(x, t), ^{x, t). The vacuum at o"_ is metastable and will de- 
cay via tunneling, and the ground state energy is E = V{a^) = O.One 
then solves a stationary functional Schrodinger equation 

S_ __.^ 



Hip = 0; where tTo- = —i-^, t^^ = —i- 



Assume the WKB ansatz ip = ^-^iplip] where 5*0 is the action of the 

bounce a (following [7] Vilenkin) and the solution for a in the thin wall 
approximation is: ( [1,7,8,10], etc) a{r) = —a tanhyu(r — _R). Here /i = 
avA, r = Vt^ + x"^- The thickness of the wall is L ^ -. R is the radius 
of the bubble at the moment of nucleation and ^R ^ 1, where R = 
^, and 5*1 = /"^ J2V{a)da is the action for 1-dimensional instanton. 
When solving for the a field, one reduces the problem to a 1 — dim QM 
problem by integrating out the other degrees of freedom (x-dependence) 
and ending up with: ^ d^x^—gH{a) = H{s),s-the parameter that 
plays the role of a time-coordinate. In this case (see Vashapati and 
Vilenkin [7]): 

So = -2MsU{s), Ms = Jd^x(^], U{s)= j d^'xiVaf + V. 

From the 0(4) symmetry under consideration, we have 

/ 2 , 2x / 2 ,2\ T —^ euclidean time 
a = a{x +T ) = a{x -t ), , , , . 

t ^ real time 

The parameter "s"is naturally chosen to be the variable : s^ = x'^+t'^ = 
x"^ — t^-. t is defined as^ = -ri-^:E = 0: and r is the Euchdean time 

' at Ms as ' ' 

. The tunneling rate for o" is 7 = e~ "b/' where S'o^, is the action in 
Euclidean time (and where t -^ — t). The well known result for the 
bounce action [1,19], is: Sq^ = 277r^S'f/2e^. 

2.3 Coordinates: Taking into consideration the 0(4) symmetries we will 
use the following coordinate system ([9]) 

rra = p sinh x 

tm = p cosh X 

Q = Q 

< p < 00,0 < X < 00 
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where rm,tm are the Minkowski radius and time. The Minkowski hne 
element written in the new coordinates is: ds"^ = —dp^ + p'^[dx^ + 
sinh^ xd'^r] . In the Euchdean region, under the barrier, the coordinates 
are: 

' PE -^ -ip 

ds% = dp% + p%{dx% + sin^ XEd'^r) 

oo> pE>0,XE^iX 

Then the future Lorentz timehke region in the above coordinates is 
the Milne universe and a = cr(p).[24]. This variable p is, clearly, the 
variable s of the previous Section 2.2. 

(2.4) Boundary Condition and Particle Creation in Functional 
Schrodinger and 2nd Quantization Picture. 

The field a interpolates under the barrier between the 2 values, (T_, (T+, 
and oo > pe > 0, p% = x"^ + t'^. We want a to reach false vac- 
uum when pe ^> oo, T ^> — oo, and the fluctuation field ip to die 
out there (i.e. at p^; ^ oo the field is in the false vacuum (o"_) 
ground state), i.e. ~ e^^^+^-'^ where — oo < r < 0, or equivalently 
~ e^'^-P'^ for pe — > oo. The In and Out frequency is defined as usual 
by c<j_|_/_ = [/c^ + m^(a+/_)]"^/^.The field reaches the turning point at 
T = t = and from then on has damped oscillations around the true 
vacuum a+ outside the barrier, i.e. 

at Pe = P = 0(r = t = 0), u — ;► (J+ (where < t < oo). 

When t — !► oo, 

\ a = a+ 

Therefore we analytically continue the solution for ip at the turning 
point Pe = P = 0, and compare it with what should have been a 
positive frequency wave in Minkowski region around a ^ cr+, i.e. with 

^ _ ^ikx-u^+t (^2.2) 

In general the analytically continued solution differs from (2.2) in that 
it contains a mixing of positive and negative frequency. This is to 
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be expected since the analytically continued solution from under the 
barrier contains the dynamics of the cr-field interpolating between 2 
values and thus the phase change is reflected and carried over when one 
continues that solution outside, in Minkowski region. These excitations 
are interpreted as particle creation. 

In the second quantization picture, the existence of a mass-squared 
term that is "time(p)-dependent" would produce the phenomena of 
particle production. Thus (^ is in a" time-dep" . background of a and 
it has been shown by [7,9] that both, functional Shrodinger and 2nd 
quantization are equivalent methods in the calculation of particle cre- 
ation. 

3. Solution for the Fluctuation Field <l>(p, %, Q) 

In the Euclidean region,the equation for the fluctuation fieldQ, which is de- 
rived from the hamiltonian in Section 2.2 is 

Pe"'^^ Pi 
where 

Li = — 5 d^ fsin^ xA ) H :^ 

In the Minkowski region the operator is —L^ = . h d^ f sinh x^x) + .■ S ■ 
Also in the Euclidean region p takes discrete values n:om to oo while in the 
Minkowski region the spectrum becomes continuous [21]. We can write 

$(p„X,,f^) = %^-F,Ux.,^) (3.2) 

PT 

where the operator L\ has eigenfunctions -Fp^m such that 
Then, after replacing (3.2) in (3.1) we obtain: 



[dl + —d,, + -T-LI- ^\<ym = (3.1) 



dl + t^JM _ 2/.2[3tanh^[/.(p, - R)] - 1] 



P 



E 



^ = (3.3) 



^Note that <I> is now used to denote the fluctuation field Lp of the previous sections. 



We have used, 



mMp,)] = m\p,) = V"[a{p^)] = 2^2[3tanh2[^(p^ - R)] - 1] 



and 



(t(p^) = — atanh[yu(pg — R)] 

as shown in Section (2.2) Equation (3.3) is a Schrodinger equation of a par- 
ticle of zero energy in a potential (E — Vq) = ^^ ^ ' ' — 'm^(p^).Thus the 
inverted potential Vq is 



Vo 



jp' + 1/4) 
Pi 



+ m^i^Pj^) (for E = Q particles). 



(3.4) 



This potential is shown in Figure 3.1. In the thin wall approximation we 
have iiR ^ 1 where p = j^, and L-thickness of the wall 



Fig 3.1 The potential V[p] for the fluctuation field of equation (3.3). 

Looking at the 4 regions we can approximate the potential (fig. 3.1) by 
(later we show that the error made in this approximation is quite negligible) 



Vn 



4/i2 

-"3 
4/i2 

4/i2-A 



as 



-k 



R + L < p^ < oo (III) 
< L {11) 



X — 



02 



P^,<P,<R-L (/) 



(3.5) 



p/pl 



< P. < 2po 



(H) 



where: 



x = (Pe - R) 



"2 = ^, «! 



^ + 2p') 



Ur, 



as 



^-Qp^ 



2/i, Ap = p2 + i/4 



The regions of interests are region //, and H. We find the solutions in 
those regions and match them. However, we included regions /// and / 
as separate, in order to indicate that the potential is nearly constant there. 
These latter regions were thus used to check the right asymptotic behaviour 
of the solutions found in the regions of interest. The turning points are: 





2^1,2 = —± J — , PQ = \h~^- (3-6) 



We will find the exact solutions in the 4 regions and match them and their 
first derivatives at the turning point. Note that if Xp ^ (n^R)'^ for large 
values of n, then the Vq defined above or the originial potential, to a very 
good accuracy, simply becomes Vq = — f +4/i^, i.e. "-|-"-term dominates the 

Pe Pe 

tanh fip and thus the distinction between the first three regions dissapears. 
In this case our result recovers that of the Kyoto group [to within a prefactor] 
[1,9] as it should since they take the mass term, m[p\ in the equation 3.4 for 
the potential, to be a constant. For sufficiently small values of p, one has 
\ < (/i-R)^, in which case the pe dependence of m^(p£;)-term should not be 
ignored since it produces interesting effects. The solutions for the 4-regions 
are: 

Region III: 

R + L < pj^ < oo, ijjjjj = A^e'^^^E ^ u = 2p 

which dies out at pe ^ oo (approaching the false vacuum from under the 
barrier) as it should from the boundary condition. 
Region II: 

R-L< Pj^<R + L, x = pj^- R 

The equation to solve is: 

[dl-L,[x-^Y + k\]i,n = Q (3.7) 

The solution to (3.7) are parabolic cylinder functions U{b,z),;U{b,—z) (or 



Debye function: U{h,z) = D_l_i i-e. 



' 2 

ilJn = A2U{b,z)+B2Uib,-z) (31 
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, Ik 1/4 0:2. 
where: 0=-^^, z = a^' (x ) 

The reason we choose U{b, ±z) is because around the turning points {a;i, X2} 
they reduce to Airy functions i.e. have the right behavior and are solutions 
there too. (If we were to use WKB we would have Airy functions by approx- 
imating Vq with a straight line around the turning points.) The asymptotic 
form for ?7(6, z) when |6| is large (in our case that is true, |6| ~ n^R^, z\^2 ~ "^ 
becomes: 



4 2/ \{2-l/ 1 Hi (() 



where: 



e = ^, t=(4|6|)i/V, r = arccoshe-eVe-l 



Thus ■?/'// is regular at the turning points with the right behavior. When 
h^ z^ (in our case h ~ 2/i^_R^, z^ ~ TWl^^' ^^ have: 



I/(t.±.)c. ^^v _^ e^"^--' (3.9) 



where 



z^ z^ 



H,«u-^-^ + ...-o, 



When 2;^ > 1 (region III), f/(6, ±2;) ~ e^""^/^ i.e. dies o ut. Thus U{b,±z) 
has the right behavior in agreement with our boundary conditions. 

Region(I): po ^ Pe, very flat. 
Solutions are: 

^jjj = Aie-'^P^ + 5ie'^^^ (3.10) 

Region(H): < pE < po 
Equation (1.3) in this region becomes: 

[C + %-V]^A = (3.11) 

Pe 
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This can be solved exactly in terms of Bessel function. The solution to (1.10) 



IS 



ipA = Pe ^±ip{'^PPE)i Cip — any modified Bessel function. 

We choose Kelvin and Neuman functions: Ki^{pE),I-ip{pE)- When Pe ^ 
and Pe -^ —ip they become Hankel functions, i.e. regular at the origin 
Pe = p = and with the right behavior when p ^ oo. 



p'/'Hd2pp) 



p—^OO 




XiLop-{ip+\m 



in the asymptotically fiat region, i.e. recover Minkowski vacuum at the true 
vacuum side of the barrier at late times. Thus: 



ijA = C2Ki^{2ppE) + CiIi^{2ppE) 

Matching at x = Xi , Xi = piE+R+^ 



(3.12) 



02 
*3 



"^ 1 ^'iiii^i) = ^'iii^i) 
Using the asymptotic form of parabolic cylinder functions U{b, ±z) (3.9) and 

after tedious algebra we get: 



2^9 



e^'^-^^AM 



A = ,/^ ^ real z, = a\'\F^ = ^a^' 



\ 2B2 = e'^-''A^{l 
where 62 = izi = l^/kA^Jos where 



A 



A2 



(Ap + 2/x2i?2 + P 



i?2 



p^R'^-Xp^p^R'^-Xp) 



XI 



i?W(Ap + 2/i2i?2 + _^ 



Xp! {p'^R^-Xp) 



(3.13) 



(3.14) 



Define q = \\/k Thus 

62 = iqxi 

when \p <^ /i^i?^ (true in our case) then A = — = \^L and 9i = 2p{R + 



02 

as 



+ A). 
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Matching at x = X2, X2 = P2E + R+ — = —A. In the same manner as 
above we have: 

f i'l{x2) = i'll{x2) 

Then (proceeding in exactly the same way) and replace {A2, -82} in favor of 
^3: 



2Ai = A3e~^^-^i{2cosh(2^2) + sinh(2^2) ■ (^^)} 
2Bi = A^e+^'^'-^'i^^^^^) sinh(2^2); ^2 - imaginary 



(3.15) 



where 63 = 2fi{R + ^ + A) 
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Matching at pE = 2po- 



^a(Po) = ^/(Po) 
V'a(Po) = ^/(Po) 



^A = pT ■ [C2K,,{2ppe) + Ci/,^(2/ip^)] ^,„»o J^[C2e"^" + Cie'-""] 




C2 



(3.16) 



The coefficient A3 is found from the unitarity condition. If i? = | ^ P, T = 
I^P, i?, T reflection and transmission coefficient in Shrodinger picture in the 
Euclidean region then R + T = 1 means | A3 P = \Ci\^ + | C2 P . 

Now that we have the solution we can analytically continue it through 
Pe = P = (turning point) in order to obtain the Bogolubov coefficient. This 
is done by replacing pE — > —ip and using the properties of Bessel functions. 
We have: ^^^^ = [C2K,pppE) + CJ,^{2ppe)]pT ■ 

In the outside region, near p — 0, that function is 



^a(p) = [a,H\l\2pp) + a2H\l\2pp)]p''\ 
But when we match iPa{.Pe — 0)(*), we use: Pe -^ —ip)- Then; 



K,^H2pp) = f e--P/2if«(2/ip) 
I.^H2pp) = ie-^/2[if«(2/ip) +iff (2/ip)] 
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and from (*) we obtain: 



02 = 1^26-^/2 



Oi 



ie"P/2[C2 + 7ri-e-"P-Ci] 



(3.17) 



r(i){2). 



When p^ 00, HlJ^^\2pp) ~ J ^^ exp[±2/iip ^ i(ip + l/2)7r/2]. Thus 

the 1 out >-state at p — > 00 in the true vacuum, iPa{.P -^ 00) is a squeezed 
state, a mixture of positive and negative frequencies: 



i}A{p -^ 00) 



026 



-^[2Mp-^-f]^^^gi[2Mp-^-f] 



(3.18) 



iFiom. relation (3.17), (3.18) (for a more rigorous treatment in 2nd quanti- 
zation, Heisenberg picture, see Kyoto, Vilenkin [7,9]) we obtain the particle 
production rate: 



\ai\ 



^p = 21" 12 (3-19) 

TT \a2\ 

4cc;2g2cos2(A2^/=^) - 7r2e-2-P-2e3 . sin2(A2v/^^) . [cj4 + ^4 

277r^; 



4(CU2 + g2)2 cOs2(A2^/=^) + sin2(A2^/^^)cj2g2 



(*) The error made in approximating the real potential to Vq, around the turning 
points is of the order (— p)^ and higher, i.e. quite negligible. Notation: uj = 

2/i, q = \yk, 62 = ^y/a^. Both ; {g. A} are real (for Ap < {pR)^). So 

sinh(2^2) -^ -isin( ^'^ ). 



?=27?\/(V + V^' + 



A2_ 



IJ.'^R*-\p) 



[see 3.14) 



A = Rj{\p + 2p^R^ 



A2 



M4R4-Ap^(M4i?4_Ap) 



To get to formula (3.19) we used (3.13)-(3.17). The calculation and ana- 
lytic continuation are straightforward. 

Relation (3.19) contains all the features found before from other authors 
(see [8,9]) i.e. the suppression for low and high energies, namely 



n. 



p ~ e "''"^ for p — small, Up ^ e 2^"^, p — large 
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but it differs greatly in terms that we considered a dynamic baclcground 
(time-depend, mass) while others have considered a constant mass-term or a 
step-function, and this dynamics of m^(p) revealed an interesting phenomena, 
resonant particle production. 

As noted before the case Xp > Afi'^R'^ where iJ,R ^ 1 corresponds to an 
almost constant mass term studied previously by Kyoto group [8] and we 
have studied the case < Ap < ^^R^ ■ n^, (ra = 0, 1, 2, 3, . . .) when the time- 
dependence of mass is important and cannot be neglected anywhere, i.e. 
when "the curvature of the bubble" becomes important. In general, from 
expression (3.19) we find that particle production is strongly suppressed by 
the factor exp[-2^2] = exp[-4/i[i?+ g + A]] since fiR > 1 and, {^, A} ~ i 
small. Also for small momentum, i.e. p small we recover the result of the 
Kyoto group in the thin wall approximate case: Up ~ e~^^^ but we have a 
different prefactor from them, we don't have (/ii?)^-facts they find, and our 
result is finite integrated over all momenta-j9. 

The most interesting feature that comes from taking into account the 
time-dependence of the mass-term, i.e. the curvature of the bubble (back- 
ground) is that for momentum "p" comparable to the curvature scale of the 
bubble - "/ii?" we get a strong enhancement of particle production, a res- 
onance effect. In this case, Xp ~ n^/x^i?^, the exponential factor becomes 

approximately zero and n^ — ./ 2, 2^2 iil-i^T \: i~z^ m \-a i ^ - In this 
case the very high momentum (p < (/i-R)^, see figure 3.3) all the sinusoidal 
functions become hyperbolic (the argument becomes imaginary) and the par- 
ticle production rate dies out very fast as expected. 

1) Thus, for momentum less than the scale of the bubble "n^" is very 
small due to the factor "e"'*^-^" . 

2) For large momentum (p ^ n^iR) one gets again a very small "np" 
suppressed by the factor "e"^?"^" . In this case there is only one region: 
Vi^ = -T — 4/i^ and Up ~ e"^^'^, as previously found [8]. 

3) Only for momentum comparable to the scale of the bubble [p ~ nuR), 
{n takes only a few integer values) one gets a resonant particle produc- 
tion. Formula (3.19) contains cases (1) and (2) as limiting cases and 
those have been previously studied by Kyoto [8,9]. Our results agree 
with theirs within a prefactor. 
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The result is that particle production is suppressed in the thin wall 
case with the exception of those modes that can resonate with the 
bubble, and those modes are not suppressed but enhanced (not shown 
previously), i.e. modes with phase "pvr" ^ n ■ 6^ = n ■ 2fiR can see 
an almost transparent barrier. If the barrier was time-dependent i.e. 
/iR = fiR(t) one would expect intuitively a larger range of modes to 
resonate with the barrier. This has to be investigated yet, (using for 
e.g. [10] model) ("n" is an integer. Only a few values of ^^Up" give a 
strong resonance, because of the suppression term q-^'^'^p-'^^-^ competing 
with the resonance effect when p is large. See (figure 3.2, figure 3.3). 

Figure 3.1 Resonant Particle Creation Number for small values of the mo- 
mentum. 



Figure 3.2 Particle creation number for high values of the momentum. 



4. Discussion. 

If we did not distinguish between fluctuation and classical fields, fiuctuations 
are the ones that contribute to the phase shift of the wave functional and give 
rise to particle production at the same time. Particle production is calculated 
in real time, after the analytic continuation of DEP outside. Thus it is the 
total phase shift coming from tunneling/Dii^P that gives rise to particle 
production. Further, on support of the above, notice that the parameter "s" 
plays the role of time. Under the barrier s < 0, outside s^ > where 



J £xH[aix, s)] = His) = ^s' + f/(s 



1 .2 , rr/ N I S = X^ +1" 
M, >0 
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Also, in a QM picture "s" is an operator, the coordinate of the particle 
in the potential U{s) thus the "kinetic energy"-(s^) is negative under the 
barrier. It has been shown that letting s run from s — » oo, to s ^ under 
the barrier is equivalent to summing up over all instantons in the dilute 
gas approximation. But instantons are paths in configuration space. Thus 
"time" -evolution of the field (i.e. let "s" run from oo to 0) is equivalent to 
the trajectory over the ensemble of paths in configuration space, (think of 
it as phase space). The changing of sign of s^ from s < to s^ > would 
correspond to an inflection point in this trajectory. That is also the turning 
point from under the barrier to outside. Taking this statistical point of view 
one would go from Euclidean to Lorentz region at the inflection point. In 
general s = s(t) is a complicated function of time, thus t^ — ^ it does not 
always give the analytical continuation from DEP to real time, i.e. is not 
crucial for tunneling, but if one can reduce the problem to a few degrees of 
freedom (s-in our case) going from "negative kinetic energy of these degrees" 
to "positive kinetic" i.e. from s^ < to s > is inevitable in the case 
of tunneling, i.e. its intrinsic to this process. The reason is because this 
parameter is defined in terms of the hamiltonian and the field is not an 
eigenstate of the hamiltonian in the case of tunneling . Outside the barrier 
s^ > and the "time" -evolution for < s < cxd gives the in-out states. 
In this approach one considers instantaneous tunneling in real time, i.e. at 
s = (in real time parameter), the field hops from the false vacuum to the 
other side of the barrier (true vacuum) and at s -^ oo settles in the true 
vacuum (1 out>-state). The solution at the turning point s = which is the 
endpoint of DEP is analytically continued thus the phase shift from DEP. 
So, particle production comes from the phase shift and the 1 in-out> states 
are obtained from the analytic continuation of the DEP in Lorentz region, 
i.e. in terms of "s"-parameter. The in-out states are the evolution of DEP 
when s ^ oo, 0. 

Then one could study real tunneling through particle production meth- 
ods. Knowing particle production and the 1 out>-state one has all the infor- 
mation about the phase shift/i.e. tunneling, and evoluting the lout>-state 
back in "time"-s one can obtain the DEP by analytically continuing at the 
inflection point (turning point where s^ goes through zero). 

Particle production calculated in our case, when the mass is time-dependent, 
i.e. dynamic background of the bubble, exhibits a strong resonance effect 
with the bubble for those modes that are comparable to the curvature scale 
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of the bubble (i.e. for pir ~ nO^ or p ~ nfiR, n = 1, 2, 3 . . . 10 . . .). When 
the modes have a high momentum they are exponentially suppressed by a 
factor Exp[—2'7ip]. This feature of exponential suppression by Exp[—27ip] 
or Exp[—4fiR] has been previsouly found by [8,9] when they considered m^- 
term to be a step function or a constant in each region respectively. In these 
last 2 limiting cases rup has a thermal behavior as previously discussed by 
[8,9]. See figure 3.2, 3.3. 

One interesting case, currently under investigation, is when the barrier 
itself is time- dependent. From naive arguments we would expect a larger 
range of modes to resonate with the bubble (as argued by Biitiker and Lan- 
duer). However we may find a nontrivial answer, different from the above, 
once that calculation is done. 

The modes of the (/9-field, if one considers those as the environment, will 
have an effect on the a-tunneling field (considered as the system). That can 
be found through calculating the noise and dissipation kernels in the case 
when the 0(4) symmetry is broken i.e. when we have more than one degree 
of freedom for a, thus bringing correlation of the WKB branches under the 
barrier. We hope to reveal these results within the next year. 
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